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ABSTRACT 



We report the discovery of a planetary system around HD9446, performed from radial velocity measurements secured with the 
spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory during more than two years. At least two planets 
orbit this G5V, active star: HD 9446b has a minimum mass of 0.7 Mj^p and a slightly eccentric orbit with a period of 30 days, whereas 
HD 9446c has a minimum mass of l.SMj^p and a circular orbit with a period of 193 days. As for most of the known multi-planet 
systems, the HD 9446-system presents a hierarchical disposition, with a massive outer planet and a lighter inner planet. 

Key words, planetary systems ~ techniques: radial velocities - stars: individual: HD 9446 
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Among the more than 400 exoplanets known so far, most of 
them have been discovered from the reflex motion they cause to 
their host-star, which can be detected from stellar radial velocity 
wobble. Thus, accurate radial velocity measurements remain a 
particularly efficient and powerful technic for research and char- 
acterization of exoplanetary systems. They allow the statistic of 
systems to be extended by completing the minimum mass-period 
diagram of exoplanets, in particular towards lower masses and 
longer periods, as the measurement accuracy improves. 

Together with the advent of the new SOPHIE spectrograph 
at the 1.93-m telescope of Haute-Provence Observatory (OHP), 
France, the SOPHIE Consortium (Bouchy et al. 2009) started 
in late-2006 a large observational program of exoplanets search 
and characterization, using the radial velocity technique. In the 
present paper we announce the discovery of two exoplanets 
around HD 9446, from radial velocity measurements secured as 
part of the second sub-program of the SOPHIE Consortium. This 
sub-program is a giant-planet survey on a volume-limited sam- 
ple around 2000 FGK stars, requiring moderate accuracy, typi- 
cally in the range 5 - 10 m s ' (Bouchy et al. |2009l l. Its goal is to 
improve the statistics on the exoplanet parameters and their host- 
ing stars by increasing the number of known Jupiter-mass plan- 
ets, as well as offer a chance to find new transiting giant planets 
in front of bright stars. SOPHIE sub-program-2 data were al- 
ready used to report the detection of several planets (Da Silva 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, by the SOPHIE Consortium (program 07A.PNPCONS). 



et al.pOOSl Santos et al.l20QS Bouchy et al. 12009] ) and to study 
stellar activity (Boisse et al. 2009^. This sub-program also aims 
at following up transiting giant exoplanets; this allowed spectro- 
scopic transits to be observed (Loeillet et al. 2008), including 
the detection of the two first cases of spin-orbit misalignment, 
nam ely XO -3b (Hebr ard et al. [2008]l then HD 80606b (Moutou 
et al. 120091 Pont et al. 120091 ). simultaneously with the discovery 
of the transiting nature of the planet in this last case. 

The SOPHIE observations of HD 9446 that allow the detec- 
tion of two new planets are presented in section|2l We derive and 
discuss the stellar and planetary properties in sections[3]and|4]re- 
spectively, and conclude in section |5] 



2. Observations 

We observed HD9446 with flie OHP 1.93-m telescope and 
SOPHIE, which is a cross-dispersed, environmentally stabilized 
echelle spectrograph dedicated to high-precision radial veloc- 
ity measurements (Perruchot et al. 120081 Bouchy et al. 2009])- 
Observations were secured in high-resolution mode, allowing 
the resolution power A/AA = 75000 to be reached. The spec- 
tra were obtained in three seasons, from November 2006 to 
March 2009. Depending of variable atmospheric conditions, the 
exposure times ranged between 3 and 18 minutes and signal-to- 
noise ratios per pixel at 550 nm between 32 and 94, with typical 
values of 5.5 minutes and 55, respectively. Exposure time and 
signal-to-noise ratio were slightly greater during the first season 
of observation. Three exposures performed through too cloudy 
conditions were excluded from the final dataset, that includes 79 
spectra. The total exposure time is about 7 hours. 
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The spectrograph is fed by two optical fibers, the first one be- 
ing used for starlight. During the first season the second SOPHIE 
entrance fiber was fed by a thorium lamp for simultaneous wave- 
length calibration. Thereafter we estimated that wavelength cal- 
ibration performed with a ~ 2-hour frequency each night (al- 
lowing interpolation for the time of the exposure) was sufficient, 
and that the instrument was stable enough to avoid simultaneous 
calibration for this moderate-accuracy program. So for the sec- 
ond and third seasons, no simultaneous thorium calibration were 
performed, avoiding pollution of the first-entrance spectrum by 
the calibration light. The second entrance fiber was instead put 
on the sky; this allowed us to check that none of the spectra 
were significantly affected by sky background pollution, espe- 
cially due to moonlight. 

We used the SOPHIE pipeline (Bouchy et al. |2009l l to exti-act 
the spectra from the detector images, cross-correlate them with 
a G2-type numerical mask, then fit the cross-correlation func- 
tions (CCFs) by Gaussians to get the radial velocities (Baranne et 
al. II 9961 Pepe et al. l2002l l. Each spectrum produces a clear CCF, 
with a 8.47 ± 0.05 km s"' full width at half maximum and a con- 
trast representing 38.0±0.8 % of the continuum. Only the 33 first 
spectral orders of the 39 available ones were used for the cross 
correlation; this allows the full dataset to be reduced with the 
same procedure, as the last red orders of the HD 9446 spectra 
obtained during the first season were polluted by argon spectral 
lines of the simultaneous wavelength-calibration. 

The derived radial velocities are reported in Table [1] The 
accuracies are between 5.4 and 9.7 ms"', typically around 
6.5ms"'. This includes photon noise (typically ~ 3.5 ms"'), 
wavelength calibration (~ 3 ms"'), and guiding errors (~ 
4 ms"') that produce motions of the input image within the fiber 
(Boisse et al. |2009bl ). These computed uncertainties do not in- 
clude any "jitter" due to stellar activity (see below). 

Table 1. Radial velocities of HD9446 measured with SOPHIE 
(full table available electronically). 



BJD 


RV 


±10- 


-2400000 


(kms"') 


(km s"') 


54043.4849 


21.7152 


0.0063 


54045.4850 


21.7451 


0.0055 


54047.4598 


21.7747 


0.0056 


54048.5033 


21.7969 


0.0054 


54889.2599 


21.7201 


0.0073 


54890.2628 


21.7409 


0.0064 


54893.2986 


21.7222 


0.0090 


54894.2961 


21.7074 


0.0097 



3. Stellar properties of HDg446 

We used the 50 SOPHIE spectra secured without simultaneous 
thorium exposure to obtain an averaged spectrum, and we man- 
aged a spectral analysis from it. Table |2] summarizes the stel- 
lar parameters. According to the SIMBAD database, HD 9446 
(HIP 7245, BDh-28 253) is a V = 8.35, high proper-motion G5V 
star. Its Hipparcos parallax {n - 19.92 + 1.06 mas) implies a dis- 
tance of53±3p c. The Hipparcos color is B-V = 0.680 + 0.015 
(Perryman et al. 11997!). 

From spectral analysis of the SOPHIE data using the method 
presented in Santos et al. ( 120041 ). we derived the temperature 



Table 2. Adopted stellar parameters for HD 9446. 



Parameters 


Values 


m,, 


8.35 


Spectral type 


G5V 


B-V 


0.680 ±0.015 


Parallax [mas] 


19.92 ± 1.06 


Distance [pc] 


53 ±3 


I' sin;* [kms"'] 


4± 1 




-4.5 ±0.1 


[Fe/H] 


0.09 ± 0.05 


Teff [K] 


5793 ± 22 


log^ [cgi] 


4.53 ±0.16 


Mass [Mg] 


1.0 ±0.1 


Radius [Rq] 


1.0 


Luminosity [Lq] 


1.1 



Teff = 5793 + 22 K, the gravity \ogg = 4.53 ± 0.16, [Fe/H] = 
+0.09 + 0.05, and = 1.0 ± 0.1 Mq. The 10 % uncertainty on 
the stellar mass is an estimation, systematic effects being difficult 
to quantify (Fernandes & Santos |2004| |. We derive a projected 
rotational velocity vsin/* =4+1 kms"' from the parameters 
of the CCF using the calibration of Boisse et al. (in preparation), 
which is similar to that presented by Santos et al. (2002). We 
also obtained [Fe/H] = +0.12+0.10 from the CCF, which agrees 
with, but is less accurate than the metallicity obtained from our 
spectral analysis. 
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Fig. 1. H and K Can lines of HD9446 on the averaged 
SOPHIE spectra. Chromospheric emissions are detected, yield- 
ing a log 7?;^^ = -4.5 + 0. 1 . 



The cores of the large Can absorption lines of HD9446 
show small emissions (Fig. [U, which are the signature of an 
active chromosphere. Such stellar activity would imply a sig- 
nificant "jitter" on the stellar radial velocity measurement. The 
level of the Ca n emission corresponds to log R'^^ - -4.5 with 
a ±0.1 dispersion according to the SOPHIE calibration (Boisse 
et al. in preparation). For a G-type star with this level of activ- 
ity, Santos et al. (2000) predict a dispersion of the order of 10 to 
20 ms"' for the stellar jitter. According to Noyes et al. ( 119841 ) 
and Mamajek & Hillenbrand (i2008i ). this level of activity im- 
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plies a stellar rotation period f rot - 10 days. This agrees with 
our vsin/* measurement, which translates into fRot < 17 days 
(Bouchy et al. ,200 5 ). depending on the unknown inclination i^, 
of the stellar rotation axes. 



4. A planetary system around HD9446 

The SOPHIE radial velocities of HD 9446 ai-e plotted in Fig. |2] 
Spanning more than two years, they show clear variations of the 
order of 200 ms"', implying a dispersion crgy - 58 ms"'; this 
is well over the expected stellar jitter due to chromospheric ac- 
tivity (10 to 20 ms"', see above). In addition, the bisectors of 
the CCF are stable (Fig. [3] upper panel), showing dispersion of 
the order of ctbis - 20 ms"', well below that of the radial ve- 
locities. An anticorrelation between the bisector and the radial 
velocity is usually the signature of radial velocity variations in- 
duced by stellar activity (see, e.g., Queloz et al. 120011 Boisse 
et al. I2009ab . The bisectors are flat by comparison with the ra- 
dial velocities, which suggests that the radial velocity variations 
are mainly due to Doppler shifts of the stellar lines rather than 
stellar profile variations. This leads to conclude that reflex mo- 
tion due to companion(s) are the likely cause of the stellar radial 
velocity variations. 
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Fig. 2. Top: Radial velocity SOPHIE measurements of HD 9446 
as a function of time, and Keplerian fit with two planets. The or- 
bital parameters corresponding to this fit are reported in Table [3] 
Bottom: Residuals of the fit with 1-cr error bars. 



These facts were known in late 2007, after two seasons of 
SOPHIE observations of HD 9446. A search of Keplerian fits 
then produced a solution with two Jupiter-like planets, on orbits 
of 30 and 190 days of period, with low eccentricities. This solu- 
tion was thereafter confirmed by the third season of observation; 
together with the "flat" bisectors, this provides a strong support 
to the two-planet interpretation of the radial velocity variations. 

Figs.|2]and|4]show the final fit of the 851-day span SOPHIE 
radial velocities of HD 9446. This Keplerian model includes two 
planets without mutual interactions, which are negligible in this 
case (see Sect.|5]). All the parameters are free to vary during the 
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Fig. 3. Bisector span as a function of the radial velocity (top) and 
the radial velocity residuals after the 2-planet fit (bottom). For 
clarity, error bars are not plotted in the bottom panel. The ranges 
have the same extents in the x- and y-axes on both panels. 



fit. The derived orbital parameters are reported in Table [3] to- 
gether with error bars, which were computed fmmx^ variations 
and Monte Carlo experiments. 

The inner planet, HD 9446b, produces radial velocity varia- 
tions with a semi-amplitude K = 46.6±3.0 ms"', corresponding 
to a planet with a minimum mass Mp sin / = 0.70 + 0.06 Mjup 
(assuming = 1.0±0.1 for the host star). Its orbit has a pe- 
riod of 30.052 ±0.027 days, and is significantly non-circular (e = 
0.20 ±0.06). This period is longer than the stellar rotation period, 
as determined above from the log R'^^ and the v sin . A Piot- 
value of 30 days would corresponds to v sin i^, < 2 km s"' , which 
is incompatible with our data. The outer planet, HD 9446c, 
yields a semi-amplitude K = 63.9 ± 4.3 ms"', corresponding 
to a planet with a projected mass Mp sin / = 1.82 ± 0.17 Mjup. 
The orbital period is 192.9 ± 0.9 days. This is about half a Earth- 
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Table 3. Fitted orbits and planetary parameters for the HD 9446 
system, with l-cr error bars. 
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Fig. 4. Phase-folded radial velocity curves for HD 9446b (P - 
30 d, top) and HD 9446c (P - 193 d, bottom) after removing 
the effect of the other planet. The SOPHIE radial velocity mea- 
surements are presented with l-cr error bars, and the Keplerian 
fits are the solid lines. Orbital parameters corresponding to the 
fits are reported in Table|3] The colors indicate the measurement 
dates. 



year, which made difficult a good phase coverage for the obser- 
vations. As seen in the lower panel of Fig. [H the rise of the ra- 
dial velocity due to HD 9446c lacks of measurements, for orbital 
phases between 0.0 and 0.3. This implies significant uncertain- 
ties on the shape of the orbit. Circularity can not be excluded 
(e = 0.06 + 0.06); furthermore, if the orbit actually is eccentric, 
there are nearly no constraints with the present dataset on the 
orientation of the ellipse with respect to the line of sight. The re- 
sulting error bars on the longitude a> of the periastron and on the 
time Tq at periastron are thus large; there are however correlated, 
and the timing of a possible transit for this planet is better con- 
strained than Tq in Table |3] Our estimations of v sin and Prot 
allow the constraint i^, > 30° to be put. So, if we assume a spin- 
orbit alignment for the HD 9446-system, = / and sin / > 0.5; 



Parameters 


HD 9446b 


HD 9446c 


P [days] 


30.052 ± 0.027 


192.9 ±0.9 


e 


0.20 ± 0.06 


0.06 ± 0.06 




-145 ±30 


-260 ± 130 


K [ms-'] 


46.6 ± 3.0 


63.9 ±4.3 


To (periastron) [BJD] 


2454 854.4 ± 2.0 


2454510 ± 70 


Mp sin / [Mj„p] 


0.70 ± 0.06+ 


1.82 ±0.17+ 


a [AU] 


0.189 ±0.006' 


0.654 ± 0.022t 


V, [Icms-'] 


21.715 ± 


0.005 


N 


79 




reduced 


2.6 




o-Q-c [ms"'] 


15.1 




Typical RV accuracy [ms"'] 


6.5 




span [days] 


851 





t: using M,^ = 1.0±0.1Mo 



this implies projected masses that translate into actual masses 
clearly in the planetary range. 

The reduced of the Keplerian fit is 2.6, and the standard 
deviation of the residuals is ctq-c = 15-1 ms"'. This is bet- 
ter than the 58-ms"' dispersion of the original radial velocities, 
but this remains higher than the 6.5-ms"' typical error bars on 
the individual measurements, suggesting an additional noise of 
~ 13.5 ms Such dispersion is precisely in the range of the 10 
to 20 m s"' expected jitter for a G-type star with this level of ac- 
tivity (Sect. [3]). Stellar activity is thus likely to be the main cause 
of the remaining dispersion, as well as the ~ 20-ms ' dispersion 
of the bisectors. The residuals of the fits do not show significant 
anticorrelation with the bisectors (Fig.[3j lower panel), as it could 
be expected in such cases (see, e.g., Melo et al. 2007' Boisse et 
al.|2009a); this is however at the limit of detection according to 
the error bars. A few bisectors values are larger than the other 
ones. They could be due to a particularly active phase of the 
star, as they are localized in a short time interval (between late 
January and early February 2007). Excluding these outliers from 
the analysis does not significantly change the results. Finally, as 
it can be seen on lower panel of Fig. |2l the residuals are signifi- 
cantly less scattered during the first season than during the third 
one. This can be mainly explained by the higher signal-to-noise 
ratio reached with longer exposure times during the first season, 
as well as the simultaneous thorium calibration secured for the 
first measurements. 

Figure|5]shows Lomb-Scargle periodograms of the radial ve- 
locity measurements of HD 9446 in four different cases: without 
any planet removed, with one or the other planet removed, and 
with both planets removed. A similar study was performed in 
the case of BD -08°2823, another star with two detected planets 
(Hebrard et al. I2009I I. In the upper panel of Fig. |5] that presents 
the periodogram of the raw radial velocity measurements of 
HD9446, periodic signals at ~ 30 days and ~ 195 days are 
clearly detected with peaks at those periods, corresponding to 
the two planets reported above, with the same amplitudes. The 
peak at ~ 1 day corresponds to the aliases of all the detected 
signals, as the sampling is biased towards "one point per night". 
A fourth, weaker peak is detected at ~ 13.3 days. A Keplerian 
fit of this signal would provide a semi-amplitude K ^ 11ms"', 
corresponding to a projected mass of 40 Earth masses. We do 
not conclude, however, that we detect a third, low-mass planet 
within the current data. 

Indeed, firstly this 13.3-day period is near the stellar rotation 
period (~ 10 days. Sect. [3]), so it could be at least partially due 
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to stellar rotation. However, no significant peaks are detected at 
this period (nor at 30 days or 193 days) on the bisectors peri- 
odograms. We interpret this 13. 3-d signal as being more likely 
due to aliases. In order to validate this, we constructed a fake 
radial velocity dataset with the same time sampling as our actual 
data, and that includes only the Keplerian model of the two plan- 
ets found above. The periodogram of this fake dataset is almost 
identical as the one plotted in the upper panel Fig.|5] it includes 
of course the two peaks corresponding to the periods of the two 
planets, but also the peak at 13.3 days. 

In addition to the one-day peak, the window function of our 
data shows a peak at ~ 14.3 days, indicating that this inter- 
val is favored in our time sampling. The 13. 3-d signal could 
thus be mainly due to the 14.3-day alias of the 192.9-day signal 
(1/13.3 ^ 1/14.3 + 1/192.9). On the second panel of Fig.|5]is 
plotted the periodogram of the residuals after subtraction of a fit 
including the 30-day-period planet only. The peak at 193 days 
is visible, as well as these aliases at 1, 13.3, and 15.4 days 
(1/15.4 ^ 1/14.3-1/192.9). In the same manner, the third panel 
of Fig. |5] shows the periodogram of the residuals after a fit in- 
cluding the 193 -day-period planet only. The peak at 193 days is 
no longer visible, and neither are the three aliases seen on the 
upper panel. This time the peak at 30 days is visible, together 
with these two aliases due to the 1-day favored sampling, at 0.97 
and 1.03 day. The bottom panel of Fig.|5]shows the periodogram 
of the residuals after subtraction of the Keplerian fit including 
HD 9446b and HD 9446c. There are no remaining strong peaks 
on this periodogram; even the 1-day alias disappeared, showing 
that most of the periodic signals have been removed from the 
data. The remaining peaks are below 10 m/s amplitude, showing 
that the main part of the detected periodical signals in our data 
are due to the two planets. The remaining signal in the residuals 
are at the limit of detection according to our accuracy. As in ad- 
dition the stellar rotation period is close to an alias of the signal 
of HD 9446c, this makes tough any characterization of radial- 
velocity signal due to stellar activity, which is expected mainly 
at the stellar rotation period. As the stellar jitter on the radial 
velocities is of the order of 10 ms this effect on the derived 
parameters of the two detected planets is negligible. 

The residuals of the measurements secured during the second 
observational season are preferentially negative (Fig. |2l lower 
panel). This may suggest a possible additional component, with 
an orbital period of the same order or larger than the time span 
of our dataset (2.3 years). Such an additional planet could not be 
established with the available data. For a ~ 2-yr period, the pro- 
jected mass of such an hypothetic planet should be lower than 
one Jupiter mass. On the other hand on short periods, a hot- 
Jupiter is excluded in this system, the accuracy of our dataset 
being good enough to detect it if there were any. The data al- 
low planets with masses larger than 0.3 Mjup and orbital periods 
shorter than 10 days to be excluded in the HD 9446 system. 



5. Discussion 

The data we presented allow us to conclude there is a planetary 
system around HD 9446, with at least two Jupiter-like planets, on 
30 and 193-day orbits. HD 9446b has a projected mass slightly 
lower than Jupiter; it is on a 0.2-eccentricity orbit, showing that 
tidal effects were not strong enough to circularize it. HD 9446c 
is at least 1 .8 times more massive than Jupiter, and is on a nearly- 
circular orbit. The host-star of this system is slightly more metal- 
lic than the Sun, in agreement with the tendency found for stars 
harboring Jupiter-mass planets (see, e.g., Santos et al. l2005l) . 




1 10 100 1000 

period (days) 



Fig. 5. Lomb-Scargle periodograms of the SOPHIE radial veloc- 
ities. The upper panel shows the periodogram computed on the 
initial radial velocities, without any fit removed. The second and 
third panels show the periodograms computed on the residuals 
of the fits including HD 9446b or HD 9446c only, respectively. 
The bottom panel shows the periodogram after the subtraction 
of the 2-planet fit. The two vertical dotted lines show the periods 
of the two planets. 



The mutual gravitational interactions between HD 9446b and 
HD 9446c are weak. The inner planet is stabilized on its orbit by 
the strong gravity of the star. Following Correia et al. (2005), a 
simulation of the two orbits from the current solution was run 
for 10^ years, in order to estimate their evolution from mutual 
interactions. This shows no significant changes in the eccentrici- 
ties, which remain in the ranges [0.18-0.23] and [0.03-0.075], 
for HD 9446b and HD 9446c respectively. Therefore this system 
is stable for 10^ years, and it seems to be stable for longer time 
scales also. We estimated the order of magnitude of the poten- 
tial transit timing variations due to those weak mutual interac- 
tions, if any of the planets of the system does transit. For that 
purpose we performed another 3-body simulation of the sys- 
tem, assuming the masses of the planets are equal to the min- 
imum masses and that the orbits are coplanar. We employed the 
Burlisch-Stoer algorithm implemented in the Mercury6 package 
(Chambers 1 19991 1 and integrated the system for 2000 days -i.e. 
around 10 orbits of the exterior planet. We found that the interac- 
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tion between the planets produces variations in the central time 
of transits with small amplitude, that does not exceed 0.4 second 
for any of the two bodies. 

No photometric search for transits have been managed for 
HD 9446; depending on the unknown inclination / of the orbit, 
the transit probability for HD 9446b and HD 9446c are about 
2 % and 1 %, respectively. There are more than 200 exoplan- 
ets detected from radial velocity surveys with orbital periods 
longer than 50 days, so with transit probabilities of the level of 
the percent. Only one is known to transit, namely HD 80606b 
(Moutou et al. 120091 ). It is likely that at least one or two more 
of these known long-period exoplanets are actually transiting as 
seen from the Earth. Their search is challenging, as the times of 
the possible transits are not known accurately, especially a few 
years after the securement of the radial velocity data. 

Among the more than 400 exoplanets discovered so far, al- 
most 25 % are located in the ~ 40 known multiple-planet sys- 
tems. Most of them have been detected from radial velocity mea- 
surements. Additional planetary companions around HD9446 
can not be detected with the available data besides the two 
planets reported, but they are of course possible, as multiple- 
planet systems are common. For example HD 155358 (Cochran 
et al. 120071 ) has a Jupiter-mass planet on an orbit similar to 
HD 9446c, and another planet with a 530-day orbital period, or 
HD 69830 (Lovis et al. 2006 ) have two Neptune-mass planets on 
orbits similar to those of the two detected planets of HD 9446, 
and a third one on a 8.7-day orbit. More data are needed, and 
the monitoring of HD 9446 should thus be maintained. As low- 
mass planets are preferentially found in multiple planetary sys- 
tems (see, e .g., McArthur et al. 2004, Pepe et al. 120071 Mayor 
et al. 120091 ). HD 9446 should be considered for high-precision 
radial-velocity programs, despite its activity level. 

The HD9446 system presents a hierarchical disposition, 
with the inner planet being the less massive one and the outer 
planet being the more massive. Fig. |6] displays the mass - 
semi-major axis relation for known multi-planetary systems. 
The data are taken from the compilation of the Extrasolar 
Planets Encyclopedi^. Most of the known multi-planetary sys- 
tems show such hierarchical disposition, as roughly the Solar 
System. The fitted relation between those two parameters pro- 
vides a planet mass proportional to a"^ (oc a'^ for the Solar 
System). The plot suggests the slope could be deeper for sys- 
tems including low-mass planets. A positive slope could be due 
to a higher migration efficiency for low-mass planets, and/or to 
the fact that giant planets are preferentially formed at larger dis- 
tances of their host stars than low-mass planets. However, ob- 
servational biases are important here, as low-mass planets are 
easier to detected at short orbital periods from radial velocity 
variations. The semi-amplitude of the reflex motion of a star due 
to a planetary companion is proportional to ^/a x Mp sin /, so 
one could expect a a^-dependance in Fig. |6l As the averaged 
slope is lower, this could suggest there is actually no strong de- 
pendance on average between those two parameters for multi- 
planet systems. 

Fig. |6] also shows that only a few multi-planetary sys- 
tems include close-in giant planets. This agrees with Wright et 
al. (2009 ), who reported that single-planet systems show a pileup 
at 3-day period and a jump at a ^ 1 AU, while multi -planet sys- 
tems show a more uniform distribution. Still, the close-in planets 
in known multi-planet system are mainly low-mass planets. Hot 
Jupiters appear to be sparse in multi-planet systems, showing 
here again a distribution which is different from single-planet 

' http://exoplanet.eu 
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Fig. 6. Semi-major axes as a function of the projected masses for 
planets in multi-planet systems. 38 known extrasolar systems are 
plotted, with the planets of a given system that are linked by a 
solid line. The HD 9446-system is shown with filled diamonds. 
The fitted relation is plotted in dotted line (mass proportional to 
fl"'^); the dashed line show the relation. The 8 planets of the 
Solar System are also plotted for comparison (circles). Systems 
with two, three and more planets are in black, blue and red, re- 
spectively. 



systems. Only five hot- Jupiters are known to be in planetary sys- 
tems, namely HIP 14810b, ups Andb, HAT-P-13b, HD 187123b, 
and HD 217107b. Single- and multi-planet systems thus appear 
to have significant differences in some of their properties. Such 
differences may provide clues allowing a better understanding 
of the formation and evolution of those systems. Improving the 
statistic of extra-solar planets should be continued, in particular 
in multi-planet systems and with radial velocity surveys. 
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